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Spin splittings in III-V materials and heterostructures are of interest because of potential appli- 
cations, mainly in spintronic devices. A necessary condition for the existence of these spin splittings 
is the absence of inversion symmetry. In bulk zincblende materials the inversion symmetry is bro- 
ken, giving rise to a small spin splitting. The much larger spin splitting observed in quantum wells 
is normally attributed to the asymmetry of the confining potential and explained on the basis of 
the Rashba effect. For symmetrically confined wells, where the only source of asymmetry is that 
of the underlying crystal potential, the confining potential strongly enhances the spin splittings. 
This enhancement does not require the asymmetry of the confining potential but depends on the 
interplay between the confinement and the crystal potential. In this situation the behavior of the 
spin splittings is consistent with the Dresselhaus contribution. 

In asymmetrically confined wells both Dresselhaus and Rashba terms contribute. 

We present a general theory of the spin splittings of these structures based on the group theory 
of diamond and zincblende heterostructures. 



I. INTRODUCTION 



II. THE COMPUTATIONAL MODELS 



Spin is one of the most intriguing properties of sub- 
atomic particles and its explanation is among the most 
significant achivcments of quantum theory. Although the 
field of electronics has taken full advantage of electronic 
charge, spin has been relatively unexploited in any prac- 
tical application with the noteable exception of magnetic 
read heads based on giant magnetoresistance. However 
recent developments suggest that this is about to change. 
The new field of spintronicsi is steadily growing with the 
aim of taking full advantage of spin as well as charge. 
Some devices have been already proposed^*^ and some 
experimental work^ is already taking advantage of some 
basic spin properties of heterostructures. 

However the spin properties of heterostructures of IV 
and III-V materials are still not completely understood. 
In spite of recent advances;^ no detailed atomistic sim- 
ulations are yet available (to our knowledge) on the spin 
properties of the materials that form the backbone of the 
semiconductor industry. Most of the work reported in the 
literature relates to a particular form of the fc-p modeli^ 
It is only recently that this model has been able to incor- 
porate the full symmetries of these heterostructuresi^ 

We shall present a method of predicting the spin split- 
tings of these structures and we shall present results 
based on the Empirical Pseudopotential Layer Method 
(EPLM) that corroborate our claims. The often dis- 
cussed contributions of bulk inversion asymmetry (BIA) 
and structural inversion asymmetry (SIA) will be pre- 
sented as consequences of the symmetries of the het- 
erostructure and we shall show clearly that any detailed 
calculation must include both. 

Section ^] will contain a summary description of the 
computational methods. Section IIIII will focus on the 
consequences of symmetry for the spin splittings of these 
structures. The following sections contain a discussion of 
our results followed by section lvTll where conclusions will 
be drawn. 



The EPLM has already been discussed elsewherai so 
we shall only describe it here briefly, mainly to point out 
the particulars of our implementation. 

In a conventional band structure calculation a set of 
eigenenergies is calculated for a particular value of the 
wavevector, k. By contrast a scattering approach, such 
as EPLM, works with a fixed energy and /c|| parallel to 
an interface and calculates a set of solutions for k± per- 
pendicular to that interface. These solutions include ex- 
amples with real and complex k±, often called the com- 
plex band structure. Solutions for different layers may 
be combined using appropriate matching conditions to 
generate solutions for more complicated combinations of 
layers. Typically the result may be expressed as a trans- 
mission coeficient for the multi-layer system. Eigenstates 
of the system then manifest themselves as resonances in 
the transmission coefficient or as bound states decaying 
into the gap of the embedding material. 

We also consider a much simpler model in which the 
system is embedded in an infinite well. 



A. The EPM and the matching conditions 

The first step of the method is to compute the complex 
band structure in each layer of the heterostructure using 
the Schrodinger equation in fc-space: 



2m 

+V{G - G')Ss.s' 



E 



jS.S 

G,G'' 



(1) 



where represents the spin-orbit term. 



The crystal potential is treated as a local pseudo- 
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potential^ described in terms of atomic form factors: 

1/(G) = 5]«„(G)5„(G), (2) 

a 

where Sa (G ) represents the structure factor and Va [G ) 
the form factors for atom species a. 

The spin-orbit term is included in its usual formula- 
tion&iaii 
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where a represents the usual vector of Pauli matrices and 
Act the spin-orbit form factors. 

In the case of hetcrostructures the symmetry is re- 
tained in the plane perpendicular to the growth direction 
and the component of the wave vector k parallel to this 
plane is still a good quantum number, k may be decom- 
posed into k = (fc|| , ki_). The Hamiltonian can be written 
as a polynomial expansion in k±: 
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where we have used the definitions: 
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It can be showniSii^ that this equation may be recast 
as an eigenvalue problem in k±_ for fixed energy, fcy and 
growth direction. 

This eigenproblem gives all the required information, 
namely all the k± and the corresponding eigenvectors, to 
allow the wavefunction to be completely determined. 

If a complex band structure is determined for adjacent 
layers i and i + with an appropriate band offset, regular 
matching conditions can be imposed as: 
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at the interface between layers. These conditions can be 
re-expressed as matrix conditions connecting the wave 
functions in both layers. A predetermined wavefunction 
in the first layer results then in fixed coefficients across 
the structure. 



B. The infinite well models 

Using the matching conditions and the complex band 
structure information from last section it is simple to fix 
infinite well conditions at the extremities of the structure. 
Denoting by the left interface for the first layer and by 
N the right interface for the end layer we shall have: 
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These conditions in conjunction with the matching 
conditions form a set of equations whose solution is usu- 
ally expressed as a determinant^ Determining the solu- 
tion is however best tackled by singular value decompo- 
sition techniques. With this approach it is easy to de- 
termine all the energy levels for the system by analyzing 
the behavior of the singular values. 

It should nevertheless be mentioned that the method 
still suffers from all the problems described previously^ 



C. The Empirical Pseudopotential Layer Method 

The EPLM is far more general than the infinite well 
models. The numerical problems inherent in that method 
are not present and appropriate boundary offsets and ma- 
terials can be selected. 

Using the complex band structure information and the 
matching techniques a scattering matrix approach may 
be implemented 4^*^ The energy levels are determined by 
analyzing the resonances of the transmission coefficient 
across the structure. By calculating the wavefunctions at 
those energies all properties are then accessible. 

This method is extremely well suited to the study of 
general hetcrostructures as no assumptions need be made 
about its layout or its growth direction. The band offsets 
between layers are taken from experimental values. 



D. Wavefunction based calculations 

Both methods give enough information, after an initial 
energy level determination, to compute the wavefunction 
or any other observablcs. In fact for a fixed growth direc- 
tion, energy and parallel wavevector the methods supply 
a complete description of the wavefunctions. This infor- 
mation is then used to compute the relevant properties. 

In particular wc compute the parallel averaged proba- 
bility density, given by: 



(11) 



which can then be used to compute the total probability 
density: 



P = J dr±a{r±). 



(12) 
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Another useful quantity is the parallel averaged spin 
polarization given by: 



= J (fri\¥{r)<j^'i>{r), (13) 
which will then give a total spin polarization of the form: 



(14) 



The method is flexible enough to compute any other 
relevant observable if necessary. 



III. SYMMETRIES 
A. Basic definitions 

Let us start by firmly setting the scope of our work. 
We are interested in the spin physics of lattice matched 
heterostructures of diamond-like and zincblende materi- 
als. 

In the possible plethora of all these structures it is use- 
ful to separate them into categories. Firstly we consider 
the atomic layer layout. If these structures have a mirror 
symmetric atomic layer distribution we say it is a sym- 
metric structure. Examples of these are often used and 
include, for instance, layouts of GaAs in AlAs. 

We have to stress that although the atomic layer layout 
may be symmetric the atomic positions within the layers 
are such that the layers arc not strict mirror images of 
each other. 

Any layout that is not symmetric is said to be asym- 
metric. 

Another useful classification considers the sharing or 
not of a common anion in the structure. The case of 
GaAs in AlAs is a clear case of a common-anion struc- 
ture. There are however situations where this does not 
happen. For example, a heterojunction of GaSb and InAs 
is such a case. This is then said to be a no-common- 
anion structure. 

These definitions will later become important in char- 
acterizing the symmetries of the structures. 



B. Symmetries 

For bulk semiconductors the spin splittings are deter- 
mined by the symmetry of the crystal lattice. Diamond 
has point group Oh and zincblende Tj. These determine 
which terms are allowed in the Hamiltonian, whether the 
spin splittings are possible and, in that case, which form 
they have. 

In the case of heterostructures the symmetry is reduced 
and it is important to know which subgroup of the bulk 
group a particular structure has. It is however impracti- 
cal to enumerate all the point groups for every possible 



layout and every possible growth direction. We confine 
ourselves therefore to the most common cases. 

In the case of diamond-like materials the most usual 
layouts consist of layers of Si and Ge grown in the [001] 
direction. In symmetric configurations these structures 
have either point group D2d or D2h- An odd number 
of atomic layers of one embedded in the other has point 
group D2d while an even number has point group I?2h- 

Non ideal interfaces containing monatomic fluctuations 
can also produce structures with point groups C2v, Civ 
and D4\i& We shall however not consider these cases as 
the methods used in this work can only handle perfect 
interfaces. 

In the case of zincblende heterostructures with a com- 
mon anion grown in the [100] direction we conclude that 
symmetric structures have point group D2d while asym- 
metric structures have point group C2v For structures 
without a common anion we have point group C2v 

A summary of these cases is given in table ^ 



Diamond 
[001] 


Symmetric 


odd number of atomic layers 


D2i 


even number of atomic layers 


D2h 


Asymmetric 




C2v 


Zincblende 
[001] 


Symmetric 


Common-anion 


D2d 


No common-anion 


C2v 


Asymmetric 




C2v 



TABLE I: Summary of point groups for heterostructures 
grown in the [001] direction. 

However, these considerations only provide us with a 
rule-of-thumb. The point group of a particular het- 
erostructure must be determined for that particular case. 
Many arrangements with only slight alterations can be 
produced which have different point groups. For other 
growth directions similar considerations apply but the 
resulting point groups will generally be different. 



C. Symmetries and the Hamiltonian 



In any particular situation we can always consider the 
Hamiltonian as an expansion in powers of fc|| about a 
high symmetry point in the 2D Brillouin zone, usually F. 
Indeed this is the basis of the popular k ■ p approxima- 
tion. The particular terms we are interested in are those 
involving spin of the form: 



7a, /3,--- X,---'^aki3 ■ --k^- 



(15) 



where Ja,f3,- 



is a case dependent constant, CTq one 

a product of 



of the Pauli matrices and the kp ■ 
components of the vector. 

Most of these terms are not allowed by symmetry and 
may be excluded. In fact, when we are interested in the 
behavior near an extrcmum such as the F point, it is 
usually sufficient to consider the first few terms. 
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In the case of diamond-hke structures with point group 
no spin-orbit terms of the form (|15|l are aUowed in 
the Hamiltonian and hence no spin splittings should be 
observed. This is easily understood as this point group 
has constituent symmetry the inversion center. 

In the case of the point group D2d this is however not 
the case. Linear terms like: 



Uyky, 



(16) 



are possible and thus linear splittings may be observed. 

This linear contribution can be understood in terms 
of the cubic termsiLiSiiS in the bulk and has hence 
been coined the Dresselhaus or bulk inversion asymmetry 
(BIA) term. 

A simple toy model can be constructed with this term 
and it is easy to calculate the spin polarization as a func- 
tion of k\\ . This dependence for a fixed magnitude of 
k\\ plotted at regular angular intervals will henceforth be 
called a spin diagram. For this particular case the two 
possible spin diagrams are displayed in figure ^ 





FIG. 1: Spin diagrams for the Dresselhaus contribution. Mag- 
nitude of spin scaled for clarity. 

In the case of a structure with C^v the only invariant 
that can be found is the term: 



(17) 



which again allows for linear splittings. This contribution 
was understood early by Rashbci22iSi in terms of the con- 
sequences of structural asymmetry in the material. This 
term is usually called the Rashba or structural inversion 
asymmetry (SIA) term. 

The toy model can be repeated with the Rashba term 
and its characteristic spin diagram is depicted in figure|21 





FIG. 2: Spin diagrams for the Rashba contribution. Magni- 
tude of spin scaled for clarity. 



The determination of the spin diagrams is an easy 
method to visualize the symmetries as these act as sym- 
metry signatures. 

For point group C-zv both BIA and SIA terms are al- 
lowed and the spin diagram looks like a superposition of 
both. An example is shown in figure |31 





FIG. 3: Spin diagrams for a particular case of mixed BIA and 
SIA contributions. Magnitude of spin scaled for clarity. 

In this situation the angle of the spin direction at fc|| = 
{kx,0) with the [100] direction is a good measure of the 
degree of mixing. 

The zero-field spin splittings may be representedii by 
an effective magnetic field. This will produce a spin 
Hamiltonian term of the form: 



H = ^ha ■ Bctf{k\\ ), 



(18) 



where i?off(A:|| ) will depend on the magnitude and direc- 
tion of the parallel wavevector fcy . From this Hamiltonian 
we obtain an overall spin splitting given by: 



A(fc|| 



h Br 



(19) 



From our previous discussion we can already deduce 
that Bcff will have two contributions: the bulk term 
(BIA) and a structural term (SIA). 

In the case of bulk zincblende structures this term pro- 
duces a well known contributioniiiiS for small values of 
k of the form: 

Bcff = Y [^-(^y - + ky{kl - kl)y + k,{kl - kl)z] 

(20) 

where 7 is a material dependent constant. In the case of 
our structures this has been showniii^ to simplify to the 
form: 

B^iA='^-^{klJ{-k,x + kyy), (21) 

where k^.^ is the value of the confined wavevector in the 
well. This term is exactly the one predicted by symmetry 
and will produce a spin splitting Abia, that is linear in 
fc|| and isotropic. 

In the case of structural inversion asymmetry it was 
shownSfliSi that this effective field is: 
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where k± is just the unit vector in the growth direction. 
For the particular case of the [001] growth direction this 
is then: 



B 



(23) 



The spin spUtting for this case Ar , is again hnear in fc|| 
and also isotropic. 

In the general case where both terms can co-exist a 
total spin splitting is given by: 



B 



BIA 



B^ 



(24) 



which can be expressed as: 

A(fc||) = y^AIj^ + A|j - 2ABiAAflsin(20), (25) 

where 9 is the angle between fc|| and the [100] direc- 
tion. In general this spin splitting is linear in k\\ but 
anisotropic. 



of layers, which materials and which band offsets are all 
set as input. This much freedom allows simulations with 
structures whose layout is completely artificial. These 
cases are however as important as those of naturally oc- 
curring heterostructures. If the latter give us precious 
data comparable with experimental results the former 
enable us to explore every possible dependence on the 
heterostructures' defining characteristics by carrying out 
computer experiments which would not be possible in a 
real laboratory. 

We should note that the artificial cases considered are 
not so far from physical situations that render then ab- 
surd. For example the infinite well situations are attain- 
able by using wide gap materials or even insulators to 
confine the system. Varying band offsets is also feasible 
to some extent by using alloying techniques. 



IV. DIAMOND HETEROSTRUCTURES 



D. Symmetries and the simulation methods 

Both models used in this work, the infinite well models 
and the Empirical Pseudopotential Layer Method, con- 
tain atomistic information about the structure under con- 
sideration and should hence reproduce the full symme- 
tries of the cases under study. There are however some 
practical but soluble problems. 

Firstly it should be noted that the method uses only 
integer numbers of monolayers (i.e. pairs of atomic lay- 
ers) which takes some cases out of our reach. An obvious 
example is the single layer of Si in Ge. As this restric- 
tion only comes about as a simplification in the matching 
technique it is possible to remove it if any of the cases in 
this category becomes important. 

Another problem concerns simulations with common 
atoms across different material layers. We should re- 
member that the Empirical Pseudopotential Method uses 
form factors determined for each material individually 
and hence the common atom is described by different 
potentials in different layers. This may cause a further 
reduction in the symmetry which is visible in the results 
but, as we shall see, does not invalidate them. A change 
to consistent atomic form factors would only partially 
solve this problem. Even if the same form factors are 
used for the common atom, and because the algorithm 
forces us to use integer number of monolayers, we would 
have one atomic layer of these anions at a different offset. 
However if we solve the matching at integer number of 
atomic layers this problem would also be solved. 

A few notes on the particular simulation cases chosen 
are appropriate at this stage. As a simulation method 
the Empirical Pseudopotential Method is quite robust 
and powerful because the layout of the heterostructure is 
completely arbitrary: the growth direction, the number 



For the crucial case of diamond heterostructures we 
have used two distinct simulations. First we considered 
the case of a layer of Ge in infinite walls. This structure 
should also have point group I?2h and no spin splittings 
should arise. The second case is that of a layers of Ge 
in Si with an artificial offset to produce a well for holes 
in the intermediate layer. This situation also has point 
group £)2ii and no spin splittings should be observed. We 
should point out that strain effects at the interfaces have 
been completely ignored. The two physical situations are 
represented in figure^ 



Ge 

lOML 



Ge 

lOML 



FIG. 4: Well layouts for diamond-like structures. The layout 
on the left represents 10 monolayers of Ge in infinite walls 
under the valence band while the right layout represents 10 
monolayers of Ge in Si with an artificial offset to produce a 
1 eV deep well adjusted on the valence band edges. 



Both simulations were run on extremely small energy 
grids and the results showed no spin splittings at all as 
we can see in figures and El for the two situations. In 
the second situation it should be noted that the two very 
close energy levels correspond to two twofold degenerate 
energy levels and not to one energy level displaying spin 
splitting. This was verified by a computation at the F 
point were the energy levels are close but not degenerate. 
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8.75^ 



3(ev) 



8.65- 



Some claims have been put forward^ that in principle 
we could engineer spin splittings by producing a structure 
with a point group different from £)2h- As we have estab- 
lished, this corresponds to the case of an odd number of 
atomic layers in the well. This situation is also extremely 
interesting: although bulk diamond-like structures do 
not show any spin splittings, and hence the Dresselhaus 
term cannot be present, in the case of hctcrostructures 
this term is present. It is hence possible to have spin split- 
tings originating from a Dresselhaus contribution even if 
the terms arc absent in bulk material. There is then an 
alternative way to engineer spin splittings in these struc- 
tures which does not rely on the Rashba effect. However 
this case is out of the reach of our simulation methods in 
their current form. 



FIG. 5: Energy dispersion for the well layout represented on 
the left of figure H fcy = (fc^,0) in units of (^) VB 
represents the valence band edge of Ge. 



There are nevertheless other ways. To show that spin 
splittings can indeed be achieved in structures involv- 
ing diamond-like materials we tested two situations that 
break inversion symmetry. The first is a sandwich of two 
layers, one of Si and the other of Ge, in infinite walls 
with an appropriate offset to line up the valence band 
edges. The second consists in the artificial case of Gc 
sandwiched between layers of Ge but with an artificial 
asymmetric offset. These cases do not have point group 
Z?2ii and do not have an inversion center; spin splittings 
are therefore allowed. Both situations are represented in 
figured 



3(ev) 



VB_Si 
0.03 0.04 



FIG. 7: Well layouts for diamond-like structures. The layout 
on the left represents 5 monolayers of Ge and 5 monolayers of 
Si in infinite walls with lined-up valence bands while the right 
layout represents 10 monolayers of Ge in Ge with artificial 
asymmetric offsets. 



FIG. 6: Energy dispersion for the layout represented on the 
right of figureH fcy = (fc^,0) in units of (^) A"\ VB_Ge 
represents the valence band edge of Ge and VB_Si the valence 
band edge of Si. 



These results are hence in line with our predictions. 



The resulting energy dispersions for these two situa- 
tions are in figures |H1 and |3 respectively. 
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8.75^ 



E(ev) 



0.02 0.03 



ations and hence these materials show potential for spin- 
tronic devices. 

The situation of infinite potential at the boundaries 
is also a case of perfect symmetric confinement. In this 
situation no extra asymmetry, other than the crystal po- 
tential, can cause the spin splittings and only an interplay 
between crystal potential and confinement is responsible 
for them. This case may be described as spin splitting 
enhancement by symmetric confinement and would have 
occurred even if a single zincblcnde compound had been 
used, as we shall see later. This enhancement refers to 
values of spin splitting far bigger than those observed in 
bulk. 

Two main conclusions may be drawn. Firstly, the two 
simulation methods we have used are perfectly capable 
of handling situations where by symmetry no spin split- 
ting is possible. Secondly, and more importantly, even 
in the case of materials where a center of inversion is 
present structures can be engineered which have spin 
splittings. This is of the utmost importance as Si and 
Ge arc presently the basis of most of the semiconductor 
industry. 



FIG. 8: Energy dispersion for the well layout represented on 
the left of figure 13 fcy = {k^,0) in units of (^) A"\ VB 
represents the artificially set common valence band edge. 



ZINCBLENDE COMMON-ANION 
HETEROSTRUCTURES 



9.H 



9.08 



9.06 



9.04 



E(ev) 



9.02 



?.98 



VB_Ge_2 



0.01 0.02 0.03 0.04 

k_x 

FIG. 9: Energy dispersion for the well layout represented 
on the right of figure |7| fcy = {kx,0) in units of [^) A 
VB_Ge_i represents the artificially set valence band edges for 
the i"' layer. 

As we can see spin splittings are apparent in both situ- 



The case of symmetric structures is still somewhat 
controversial in the literature. Most of this controversy 
stems from the fact that the conventional model of elec- 
tronic structure in hetcrostructurcs, the k ■ p method, 
does not fully account for their symmetries. In fact, in 
this method the wavefunctions are expanded in a set of 
F Bloch states of the zincblende crystal. Further this 
expansion is restricted to a few states, usually the top 
of the valence band and the bottom of the conduction 
band. With this set it is impossible to resolve any atom- 
istic details and the method is incapable of reproducing 
the correct point group symmetries of the structure. This 
was thought not to be problematic given the small values 
of the bulk terms from BIA. However, worki^ as early 
as 1988 hinted that this is not the case. More recent 
theoretical studies^ have confirmed this. Nevertheless 
all these studies rely on introducing terms in the k ■ p 
model that mimic the symmetries of the structure under 
consideration and have thus to be tailored to particular 
situations. In contrast any atomistic approach, like the 
Empirical Pseudopotential Method, incorporates by con- 
struction the correct symmetry of the structure. How- 
ever no such calculations, or even experimental data, is, 
to our knowledge, available for the case of spin split- 
tings. It should also be noted that in this particular 
case, zincblende common-anion structures, no structural 
asymmetry is introduced and any spin splitting cannot 
be attributed to the Rashba effect. 

Given that linear or cubic terms exist in bulk it is ex- 
pected that any band would split linearly close enough 
to the F point. The fact that the point group of these 
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structures with a common-anion D2d allows these terms 
further reinforces our belief that this must indeed hap- 
pen. Nothing however tells us that the coefficient associ- 
ated with this phenomenon would be big enough to allow 
eventual technological use of these structures. Neverthe- 
less the case of spin splitting enhancement by symmetric 
confinement that we have already encountered lets us be- 
lieve that this is the case. 

A first introductory calculation with GaSb in infinite 
walls is then performed. The band edge layout for this 
structure is depicted on the left of figure (|10|l . The en- 
ergy dispersion computed with the infinite well model is 
represented in figure ITTl 



case of symmetric confinement enhancement of the spin 
splittings. It should be noted that for this particular di- 
rection there is no spin splitting in the bulk case. 



A more realistic case of GaSb sandwiched between 
AlSb was also used. The band layout is depicted on the 
right of figure ^| where the band offsets were set to ac- 
knowledged experimental valuesiS^ 



GaSb 
lOML 



GaSb 
lOML 



0.35 eV 
AlSb 



The computed energy dispersions for both the conduc- 
tion and valence band energy windows is shown in fig- 
ures ^1 and ^1 respectively. 



FIG. 10: Well layout for zincblende structures. The layout 
on the left represents 10 monolayers of GaSb in infinite walls 
under the valence band while the right layout represents a 
well of 10 monolayers of GaSb in AlSb. 




E(ev) 



FIG. 11: Energy dispersion for the layout represented on the 
left of figureHSl fcy = (fc^,0) in units of (^)A-\ VB repre- 
sents the valence band edge of GaSb. 

Again we clearly see spin splittings that far exceed typ- 
ical values of bulk spin splittings. This is then another 



FIG. 12: Energy dispersion for the layout represented on the 
right of figure [TUl for the energy window of conduction band. 
k\\ = {kx,0) in units of [^)A~^. CB_GaSb represents the 
conduction band edge of GaSb. 
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FIG. 13: Energy dispersion for the layout represented on the 
right of figure ITI71 for the energy window of the valence band. 
fc|| = {kx,0) in units of (^)A"^ VB.GaSb represents the 
valence band edge of GaSb. 

Every band is clearly spin split: a fact that can be 
confirmed by a calculation of the spin polarization. As 
before, for the [001] direction no spin splitting is observed 
in the bulk case. In both the conduction and valence 
band the spin splitting can be easily fitted to a linear 
dispersion, AE oc k, yielding a coefficient of 76.5 meV A 
and 196.58 mcV A respectively. 





FIG. 14: Spin diagram for the first conduction level (left) 
and valence level (right) in the band structure of figures 1121 
& [SI respectively with fc,, = 0.02 (^)A~\ fcy = {kx,ky) in 
units of (^)A~^. Magnitude of spin scaled for clarity. 



this case the spurious Rashba term is approximately 0.34 
times the Dresselhaus contribution. It should be noted 
that the z component of the spin polarization is always 
found to be zero within numerical fluctuations. Spin dia- 
grams for other energy levels were computed with similar 
results. 

It is important to note the values of the linear coeffi- 
cients. These values are comparable to the linear coeffi- 
cients stated in literatures^ for the Rashba coefficient in 
asymmetric structures. It should also be noted that these 
values cannot be attributed to the spurious Rashba con- 
tributions in this case. The situation with infinite walls 
does not suffer from this contamination and produces 
similar results. This fact alone is technologically impor- 
tant: structural asymmetry is probably not required to 
produce structures that behave similarly to those cur- 
rently proposed for the purpose of creating spin split- 
tings. The spin behavior as shown in the spin diagrams 
is completely different, however. 

As the SIA and BIA contributions have similar magni- 
tudes it is important that both are included in any study. 
This is particularly important when methods, such as k-p 
are employed, which don't automatically contain Dressel- 
haus contributions. Even when the Rashba contribution 
is important the interplay of the 2 terms may produce 
sizable effects. 

Another major achievement of this method is the pos- 
sibility of extracting atomistic details in clear contrast 
to the majority of the methods previously used. In fig- 
ure El the parallel averaged probability densities for the 
first energy level in the valence band energy window is 
depicted. 




The spin diagrams for a fixed k\\ = 0.02 (^)A~^ and 
computed for the first energy level on the energy windows 
for the conduction and the valence bands are shown in 
figure^! They represent a clear signature of the Dressel- 
haus terms and cannot be attributed to any structural in- 
version asymmetry. The slight deviation from the perfect 
Z?2d signature was already explained in section Fill Dl In 



FIG. 15: Parallel averaged probability density for the first 
energy level in the energy window of the valence band, as 
depicted in figure [T^ for fcy = (0.02,0) (^)A~\ Horizontal 
axis corresponds to the growth direction in Angstrom. 

Clearly the general trend is the same as in previous cal- 
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culations but a wealth of extra information is portrayed. 
The general belief that the envelope behaves as predicted 
by "particle-in-box" type calculations is confirmed even 
in this atomistic calculation. Also displayed clearly is 
that there are deviations from the envelope behavior in 
the atomistic detail of these graphs. This could help en- 
gineer particular structures tailored to exhibit particular 
physical effects or even in the determination of the best 
doping technique. Similar results arc obtained for all 
other energy levels. 

Atomistic detail is also clear in the parallel averaged 
spin polarization. This is depicted in figure 1161 corre- 
sponding to the parallel averaged probability density de- 
picted in figure El This level of detail in spin behavior 
can be important in tailoring particular spin properties 
and possibly in doping with magnetic materials. 




FIG. 16: Three-dimensional representation of the parallel av- 
eraged spin polarization S{ri_) = {ax,o'y,a',), corresponding 
to the parallel averaged probability density in figure IT^ Ver- 
tical axis corresponds to the growth direction in Angstrom. 
The continuous line corresponds to following the tip of the 
vector (T{r±) in space after appropriate scaling. The arrow is 
a vector in the direction of fcy introduced as guidance. 



It should be noted that most of the atomistic detail, 
in the particular case of the Dresselhaus term, is in a 
plane perpendicular or nearly perpendicular to k\\. It's 

however the component parallel or nearly parallel to fc|| 
that averages to the total spin polarization in accordance 
to the spin diagrams characteristic of this term. This 
atomistic detail of the spin polarization has never been 
reported previously. 

It is also possible to use known data from bulk mate- 
rials to compare with our results. We know that for the 
conduction band spin splittings are given by cubic terms 
originating in the Hamiltonian term given by equation 
(|2l)|l . The constant 7 can be obtained in the literature24 

and values range from 109.4 to 153.9 cV A for theoreti- 
cal predictions with several methods and 186.3 eV A for 
the experimental value. In the case of GaSb in AlSb the 
spin splittings of the levels in the conduction band energy 
window should then originate, to first approximation, in 
linear and quadratic terms in fc^ the value of the con- 
fined wavevector. For wide enough wells this value should 
be, for the first energy level, approximately j^, with L the 



well width given by \aN where a is the lattice constant 
and N the number of monolayers in the well. 

A well width dependence for the first electronic en- 
ergy level was then computed and the result is shown in 
figure [T7I 



FIG. 17: Spin splitting of the first energy level in the energy 
window of the conduction band for fcy — (0.01, 0) (^) A~^ as 
a function of the inverse well width given by the number of 
monolayers of GaSb. 

For the widest wells the linear contribution should be 
dominant and the spin splitting should behave as: 



(26) 



A linear fit to this data gives a value of 7 of 118.2 eV A . 

Given all the approximations used this value is in very 
good agreement not only with results using very different 
approaches, but also with experiment, and gives us con- 
fidence in the reliability all of the results of the method. 

As a final conclusion on this section it should be noted 
that although our results are for a particular case we be- 
lieve that the dominant trends in these structures are 
determined by symmetry and rather than by the partic- 
ular atoms involved. Further calculations for the case of 
GaAs in AlAs were also computed and produced similar 
results. 



VI. ASYMMETRIC HETEROSTRUCTURES 

In the case of symmetric wells with a common anion we 
probed the consequences of introducing the Dresselhaus 
term. As a first approach we would be interested in prob- 
ing the Rashba contribution in the same way. The case 
of asymmetric structures grown in the [001] direction is 
however usually in the C2V point group class which allows 
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both terms. However wc already know that the Drcssel- 
haus term arises from bulk inversion asymmetry. Hence, 
if we consider a structure with dominant structural inver- 
sion asymmetry we should minimize the its effects. For 
this purpose we revisit the artificial structure of Ge sand- 
wiched between layers of Ge but with asymmetric band 
offsets as depicted on the right of figure [7| Wc should 
nevertheless remember that, even though bulk Ge does 
not allow the Dresselhaus contribution, this structure will 
contain this term even if it is small. We have already 
shown that such contributions occur in cases where the 
center of inversion has been removed as in the case of the 
double layer of Si and Ge in infinite walls. However these 
should be smaller than in any structure constructed out 
of zincblcndc materials. As we shall see the Dresselhaus 
contribution in this case is far smaller than the Rashba 
term. 

The computed energy dispersion in the energy window 
of the valence band was already displayed in figure 
The bands clearly show spin splittings which have been 
confirmed by a calculation of the spin polarization. A 
linear fit to the actual spin splitting for the first band 
gives a splitting coefficient of 324.3 meV A. 

This value is of the same order as those calculated for 
the case of symmetric structures; further reinforcing our 
conclusion that BIA must always be taken into consider- 
ation. 




FIG. 18: Spin diagram for the first valence level with k\\ = 
0.02 (^)A~\ = {ka:,ky) in units of (^)A"\ Magnitude 
of spin scaled for clarity. 



More importantly wc have also calculated the spin dia- 
gram for this case which is given in figure lT^ In this case 
the z-component of the spin polarization is also found to 
be zero within numerical fluctuations. As we can see 
it forms a clear signature of the Rashba contribution. 
The slight deviation is due, in this case, to the Dressel- 
haus term which is allowed by the symmetry arrange- 
ment. This term is approximately 0.14 times the domi- 
nant Rashba contribution. 

In this situation it is also possible to extract atomistic 
details from our results. In figure [T^ a parallel averaged 
probability density is shown. 




O lO 20 30 



FIG. 19: Parallel averaged probability density for the first 
energy level in the energy window of the valence band, as 
depicted in figureEl for fcy = (0.02,0) (^)A"\ Horizontal 
axis corresponds to the growth direction in Angstrom. 

The typical envelope behavior is reproduced but sig- 
nificantly more information is present. As in the pre- 
vious case this might be significant for engineering new 
structures. Atomistic detail is also present in the parallel 
averaged spin polarizations. This quantity is depicted in 
figure I2UI corresponding to the parallel averaged proba- 
bility density shown previously. 




FIG. 20: Three-dimensional representation of the parallel av- 
eraged spin polarization ff{r±) = (cto; , a ^ , ) , corresponding 
to the parallel averaged probability density in figure [191 Ver- 
tical axis corresponds to the growth direction in Angstrom. 
The continous line corresponds to following the tip of the vec- 
tor ff{r±) in space after appropriate scaling. The arrow is a 
vector in the direction of fcy introduced as guidance. 

This case is different from the case of the Dresselhaus 
term. Most of the atomistic detail is in the direction 
of fc|| while the perpendicular component averages to the 
total spin polarization in accordance to the spin diagrams 
characteristic of the Rashba term. Again this level of 
atomistic detail has never been reported previously. 
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The main conclusion of our calculations is that struc- 
tural asymmetry in the system introduces a new term in 
the Hamiltonian that allows linear splittings close to the 
r point in perfect agreement with the Rashba effect. It is 
also important to emphasis that the Empirical Pseudopo- 
tential Method is capable of simulating this situation. 

It is possible to minimize the BIA contribution by us- 
ing diamond-like materials rather than producing struc- 
tures with point group C^v However in a general case 
both terms will interplay and the magnitude of the Dres- 
selhaus contribution can be comparable to that of the 
Rashba term. A detailed treatment of these structures 
must always incorporate both to guarantee an accurate 
description of the physical processes. This was previ- 
ously discussed in the literature^ using the independent 
k ■ p method. 



VII. CONCLUSIONS 

We have shown that the phenomenon of spin split- 
tings in heterostructures can be put into a consistent 
global framework. Until now the physics of symmetric 
and asymmetric structures was thought to be in essence 
different. We have proven that they are however just dif- 
ferent expressions of the same underlying physics: "sym- 
metry rules" . 

The Empirical Pseudopotential Layer Method is very 
well suited for atomistic detailed calculations for these 
structures. It can accurately predict the energy levels 
and their spin splittings throughout the Brillouin zone. 
The case of the neighborhood of the P point, specifically 
analyzed in this chapter due to its technological impor- 
tance, is just a particular case. The method can also 
play a crucial role in determining both transport and op- 
tical properties of these structures as the wavefunction 
for every possible state is easily obtained. It is also flex- 
ible enough to handle arbitrary growth directions and 
automatically incorporating the correct symmetries, in 
contrast to the k ■ p method which requires to be adapted 
for each particular symmetry case. 

Our most significant conclusion is that in all the stud- 
ied structures linear terms emerge in the spin splittings. 
These terms play a crucial part in determining the top or 
bottom of bands near the P point and are usually the de- 
termining factor in experimental results that probe close 
to that point. 

The results involving diamond-like materials also pro- 
vide substantial conclusions. Although bulk materials 
with this structure do not exhibit spin splittings we 
proved that heterostructures with these materials may 
exhibit them. This is of crucial importance as most of 
today's semiconductor industry is based on Si and Ge. 
The technology is still not capable of producing struc- 
tures with absolute layer precision and so the case of odd 
number of atomic layers is probably unlikely to be of 
use. However this apparent technological problem can 
actually be explored to produce interface roughness that 



is responsible for lowering of symmetry; thus making spin 
splittings possible^. 

Symmetric structures with a common anion exhibit lin- 
ear splittings comparable to those determined in asym- 
metric structures but with different spin behavior. The 
new field of spintronics may find this new degree of free- 
dom technologically useful. The combination of both 
Dresselhaus and Rashba may be used to modulate a par- 
ticular spin behavior. 

Our results are in good agreement with the most recent 
theoretical studjs^ in the literature. All linear coefficients 
for both valence and conduction bands are similar even 
though the methods used are different. The conclusion 
that both BIA and SIA must always be included in every 
calculation is also drawn there. 

Comparison with experiment is more difficult, how- 
ever. Firstly there is a lack of results for the diamond-like 
structures and for the zincblende symmetric common- 
anion cases due to the relevance always given to the 
asymmetric case. Secondly, even for the case of asymmet- 
ric structures, it has emerged recently that the most com- 
mon experimental procedure, using Shubnikov-de Haas 
oscillations, might not have been properly analyzed^^. 
Neveretheless the values obtained^^ are of the same or- 
der of magnitude. 

There is also another indirect experimental result that 
can be related to our prediction. An in-plane polariza- 
tion anisotropy is observed^ in the case of structures of 
C2v point group relevant to optical considerations involv- 
ing interband transitions^^. This anisotropy reveals itself 
in the band structure of these heterostructures by a clear 
difference between the [110] and the [TlO] directions. This 
difference is visible in our computed band structures for 
the asymmetric case. For the common-anion symmet- 
ric situation no appreciable difference could be detected. 
This effect is clearly important for optical devices and 
can also be determined by the Empirical Pseudopoten- 
tial Method. 

Another point usually considered^ in the context of 
spin splittings in heterostructures is the influence of the 
main gap on the relative magnitude of both contribu- 
tions. This claim states that in narrow gap systems SIA 
effects dominate while for wide band gaps BIA is larger. 
Although we do not have enough data to compare linear 
coefficients for narrow and large band gaps our results are 
enough to conclude that even in narrow gap systems, like 
GaSb, both contributions should be taken into consider- 
ation. We consider this claim to be an oversimplification 
of the dependence of spin splittings on the parameters of 
the structure. These will depend in a non-simple way not 
only on the main gap but also on the well width, depth 
and applied flelds. Results obtained^ elsewhere seem to 
confirm this. 

The results presented here are the first, to our knowl- 
edge, atomistic simulations that show that the full spin- 
orbit interaction caused by the atomic cores is the dom- 
inant contribution for the zero-field spin splittings. The 
particular symmetry of the case under consideration de- 



13 



termines the possible behavior: Dresselhaus or Rashba. 
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